Population genomics of flat-tailed horned lizards (Phrynosoma mcallii) informs conservation and management across a fragmented Colorado Desert landscape.

Phrynosoma mcallii (flat-tailed horned lizards) is a species of conservation concern in the Colorado Desert of the United States and Mexico. We analysed ddRADseq data from 45 lizards to estimate population structure, infer phylogeny, identify migration barriers, map genetic diversity hotspots, and model demography. We identified the Colorado River as the main geographic feature contributing to population structure, with the populations west of this barrier further subdivided by the Salton Sea. Phylogenetic analysis confirms that northwestern populations are nested within southeastern populations. The best-fit demographic model indicates Pleistocene divergence across the Colorado River, with significant bidirectional gene flow, and a severe Holocene population bottleneck. These patterns suggest that management strategies should focus on maintaining genetic diversity on both sides of the Colorado River and the Salton Sea. We recommend additional lands in the United States and Mexico that should be considered for similar conservation goals as those in the Rangewide Management Strategy. We also recommend periodic rangewide genomic sampling to monitor ongoing attrition of diversity, hybridization, and changing structure due to habitat fragmentation, climate change, and other long-term impacts.

The Artefactual Branch Effect and Phylogenetic Conflict: Species Delimitation with Gene Flow in Mangrove Pit Vipers (Trimeresurus purpureomaculatus-erythrurus complex).

Mangrove pit vipers of the Trimeresurus purpureomaculatus-erythrurus complex are the only species of viper known to inhabit mangroves. Despite serving integral ecological functions in mangrove ecosystems, the evolutionary history, distribution, and species boundaries of mangrove pit vipers remain poorly understood, partly due to overlapping distributions, confusing phenotypic variations, and the lack of focused studies. Here, we present the first genomic study on mangrove pit vipers and introduce a robust hypothesis-driven species delimitation framework that considers gene flow and phylogenetic uncertainty in conjunction with a novel application of a new class of speciation-based delimitation model implemented through the program Delineate.Our results showed that gene flow produced phylogenetic conflict in our focal species and substantiates the artefactual branch effect where highly admixed populations appear as divergent non-monophyletic lineages arranged in a stepwise manner at the basal position of clades. Despite the confounding effects of gene flow, we were able to obtain unequivocal support for the recognition of a new species based on the intersection and congruence of multiple lines of evidence. This study demonstrates that an integrative hypothesis-driven approach predicated on the consideration of multiple plausible evolutionary histories, population structure/differentiation, gene flow, and the implementation of a speciation-based delimitation model can effectively delimit species in the presence of gene flow and phylogenetic conflict.

A new species of Dipsas (Serpentes, Dipsadidae) from central Panama.

A new species of Dipsas Laurenti, 1768, from Central Panama is described based on molecular analyses, hemipenial morphology, and external characters. This is the sixth species of Dipsas to be described for the country; the snake has been suspected to exist since 1977 and has not been thoroughly studied until now. Additionally, morphological comparations including scale counts are done with other species within the genus, and the current geographic distribution of Dipsastemporalis (Werner, 1909), the sister species, is updated. Finally, a key to the species of Dipsas currently known from Middle America is presented.

ResumenDescribimos una nueva especies de Dipsas Laurenti, 1768 de la región central de Panamá en base a análisis moleculares, morfología hemipenial y caracteres de morfología externa. Esta es la sexta especie del género Dipsas descrita para el país. Se sospechaba su existencia desde 1977 pero no había sido estudiada exhaustivamente hasta ahora. Adicionalmente, presentamos comparaciones morfológicas (incluyendo lepidosis) con otras especies del género y actualizamos la distribución geográfica de su especie hermana Dipsastemporalis (Werner, 1909). Finalmente, presentamos una clave para las especies de Dipsas distribuidas en Centroamérica.

A new species of pit-viper from the Ayeyarwady and Yangon regions in Myanmar (Viperidae, Trimeresurus).

In a genomic study by Chan and colleagues, pit-vipers of the Trimeresuruserythrurus-purpureomaculatus complex from the Ayeyarwady and Yangon regions in Myanmar were demonstrated to be a distinct species based on robust population genetic and species delimitation analyses. Here, we provide morphological characterizations and a formal description of those populations as a new species. The new species, Trimeresurusayeyarwadyensissp. nov., is most closely related to T.erythrurus and T.purpureomaculatus and shares morphological characteristics with both of those species. Some specimens of T.ayeyarwadyensissp. nov. have green dorsal coloration and no distinct dorsal blotches (a trait shared with T.erythrurus but not T.purpureomaculatus), while others have dark dorsal blotches (a trait shared with T.purpureomaculatus but not T.erythrurus). The distinct evolutionary trajectory of the new species, coupled with the lack of obvious morphological differentiation, represents a classic example of the cryptic nature of species commonly found in the Trimeresurus group of Asian pit-vipers and underscores the need for data-rich analyses to verify species' boundaries more broadly within this genus.

DNA barcoding of the National Museum of Natural History reptile tissue holdings raises concerns about the use of natural history collections and the responsibilities of scientists in the molecular age.

Natural history collections are essential to a wide variety of studies in biology because they maintain large collections of specimens and associated data, including genetic material (e.g., tissues) for DNA sequence data, yet they are currently under-funded and collection staff have high workloads. With the advent of aggregate databases and advances in sequencing technologies, there is an increased demand on collection staff for access to tissue samples and associated data. Scientists are rapidly developing large DNA barcode libraries, DNA sequences of specific genes for species across the tree of life, in order to document and conserve biodiversity. In doing so, mistakes are made. For instance, inconsistent taxonomic information is commonly taken from different lending institutions and deposited in data repositories, such as the Barcode of Life Database (BOLD) and GenBank, despite explicit disclaimers regarding the need for taxonomic verification by the lending institutions. Such errors can have profound effects on subsequent research based on these mis-labelled sequences in data repositories. Here, we present the production of a large DNA barcode library of reptiles from the National Museum of Natural History tissue holdings. The library contains 2,758 sequences (2,205 COI and 553 16S) from 2260 specimens (four crocodilians, 37 turtles, and 2,219 lizards, including snakes), representing 583 named species, from 52 countries. In generating this library, we noticed several common mistakes made by scientists depositing DNA barcode data in public repositories (e.g., BOLD and GenBank). Our goal is to raise awareness of these concerns and offer advice to avoid such mistakes in the future to maintain accurate DNA barcode libraries to properly document Earth's biodiversity.

Interspecific Gene Flow and Mitochondrial Genome Capture during the Radiation of Jamaican Anolis Lizards (Squamata; Iguanidae).

Gene flow and reticulation are increasingly recognized as important processes in the diversification of many taxonomic groups. With the increasing ease of collecting genomic data and the development of multispecies coalescent network approaches, such reticulations can be accounted for when inferring phylogeny and diversification. Caribbean Anolis lizards are a classic example of an adaptive radiation in which species have independently radiated on the islands of the Greater Antilles into the same ecomorph classes. Within the Jamaican radiation at least one species, Anolis opalinus, has been documented to be polyphyletic in its mitochondrial DNA, which could be the result of an ancient reticulation event or incomplete lineage sorting (ILS). Here, we generate mtDNA and genotyping-by-sequencing (GBS) data and implement gene tree, species tree, and multispecies coalescent network methods to infer the diversification of this group. Our mtDNA gene tree recovers the same relationships previously inferred for this group, which is strikingly different from the species tree inferred from our GBS data. Posterior predictive simulations suggest that our genomic data violate commonly adopted assumptions of the multispecies coalescent model (MSCM), so we use network approaches to infer phylogenetic relationships. The inferred network topology contains a reticulation event but does not explain the mtDNA polyphyly observed in this group; however, coalescent simulations suggest that the observed mtDNA topology is likely the result of past introgression. How common a signature of gene flow and reticulation is across the radiation of Anolis is unknown; however, the reticulation events that we demonstrate here may have allowed for adaptive evolution, as has been suggested in other, more recent, adaptive radiations. [Adaptive radiation; hybridization; introgression; multispecies network coalescent; posterior predictive simulation.].

Evidence of two deeply divergent co-existing mitochondrial genomes in the Tuatara reveals an extremely complex genomic organization.

Animal mitochondrial genomic polymorphism occurs as low-level mitochondrial heteroplasmy and deeply divergent co-existing molecules. The latter is rare, known only in bivalvian mollusks. Here we show two deeply divergent co-existing mt-genomes in a vertebrate through genomic sequencing of the Tuatara (Sphenodon punctatus), the sole-representative of an ancient reptilian Order. The two molecules, revealed using a combination of short-read and long-read sequencing technologies, differ by 10.4% nucleotide divergence. A single long-read covers an entire mt-molecule for both strands. Phylogenetic analyses suggest a 7-8 million-year divergence between genomes. Contrary to earlier reports, all 37 genes typical of animal mitochondria, with drastic gene rearrangements, are confirmed for both mt-genomes. Also unique to vertebrates, concerted evolution drives three near-identical putative Control Region non-coding blocks. Evidence of positive selection at sites linked to metabolically important transmembrane regions of encoded proteins suggests these two mt-genomes may confer an adaptive advantage for an unusually cold-tolerant reptile.

A new cryptic species of fringe-toed lizards from southwestern Arizona with a revised taxonomy of the Uma notata species complex (Squamata: Phrynosomatidae).

Fringe-toed lizards (Uma) are among the most specialized lizards in North America, adapted to insular windblown sand habitats in the hyper-arid southwestern deserts, with allopatric distributions, subtle morphological variation, and an unstable taxonomic history. We analyzed a morphological dataset of 40 characters for 65 specimens and a molecular dataset of 2,286 bases from three mitochondrial loci for 92 individuals and interpreted these data alongside published analyses of multi-locus genetic data with the goal of revising the taxonomy of the Uma notata (Baird 1858) species complex. We confirmed that fringe-toed lizards from the Mohawk Dunes in southwestern Arizona (U. sp.) constitute a cryptic species sister to the rest of the complex that can be diagnosed with DNA barcoding and geography, so we describe and name this species Uma thurmanae sp. nov. We also confirmed the evolutionary distinctiveness of U. inornata (Cope 1895), an endangered species endemic to Coachella Valley in southern California. We designate a lectotype for the taxon U. "rufopunctata", but we put its name in quotation marks to reflect its uncertain taxonomic status with respect to its neighboring species U. cowlesi and U. notata.

The tuatara genome reveals ancient features of amniote evolution.

The tuatara (Sphenodon punctatus)-the only living member of the reptilian order Rhynchocephalia (Sphenodontia), once widespread across Gondwana1,2-is an iconic species that is endemic to New Zealand2,3. A key link to the now-extinct stem reptiles (from which dinosaurs, modern reptiles, birds and mammals evolved), the tuatara provides key insights into the ancestral amniotes2,4. Here we analyse the genome of the tuatara, which-at approximately 5 Gb-is among the largest of the vertebrate genomes yet assembled. Our analyses of this genome, along with comparisons with other vertebrate genomes, reinforce the uniqueness of the tuatara. Phylogenetic analyses indicate that the tuatara lineage diverged from that of snakes and lizards around 250 million years ago. This lineage also shows moderate rates of molecular evolution, with instances of punctuated evolution. Our genome sequence analysis identifies expansions of proteins, non-protein-coding RNA families and repeat elements, the latter of which show an amalgam of reptilian and mammalian features. The sequencing of the tuatara genome provides a valuable resource for deep comparative analyses of tetrapods, as well as for tuatara biology and conservation. Our study also provides important insights into both the technical challenges and the cultural obligations that are associated with genome sequencing.

Publisher Correction: The tuatara genome reveals ancient features of amniote evolution.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Six additional mitochondrial genomes for North American nightsnakes (Dipsadidae: Hypsiglena) and a novel gene feature for advanced snakes.

The North American nightsnakes in the genus Hypsiglena is composed of nine named and at least two unnamed species. Here, we provide the first mt-genome of H. affinis, an additional mt-genome for H. sp. nov. 1, and four additional mt-genomes from the widespread H. jani. These mtDNA genomes were sequenced using both Illumina and Ion Torrent sequencing technologies. The resulting genomes contained the expected 13 protein coding genes, 22 tRNA genes, 2 rRNA genes, and 2 control regions typical of colubroid snakes. Two of the H. jani samples had partial tRNAIle genes upstream of CR2 which has not been previously documented in colubroid snakes. A maximum likelihood gene-tree based on these data combined with previously published sequence data recovers a well-supported phylogeny and is in concordance with previous estimates of evolutionary relationships in this group.

The discovery of Rock Geckos Cnemaspis Strauch, 1887 (Squamata: Gekkonidae) in the Tanintharyi Region, Myanmar with the description of two new species.

We report the first occurrence of the Asian Rock Gecko genus Cnemaspis Strauch, 1887 from mainland Myanmar based on a series of specimens recently collected from the Tanintharyi Region. These records come from several localities in the Tanintharyi Region and fill a significant sampling gap for the genus. Molecular phylogenetic analyses using the mitochondrial gene ND2 identify two distinct clades, the first includes Cnemaspis siamensis of the Southeast Asian Cnemaspis group and the second includes two new species belonging to the South Asian Cnemaspis kandiana group. These two species are morphologically distinct and are distinguished from all other members in the C. kandiana group by a combination of character states. The first species, Cnemaspis tanintharyi sp. nov. occurs on the mainland in southern Tanintharyi and is distinguished from all congeners by the possession of keeled pectoral scales; smooth ventral scales and abdominal scales; 2-4 precloacal pores; 4-5 femoral pores on each leg; smooth subcaudal scales with the median row enlarged; coloration of the gular region beige, dark gray-brown with dark blue hueing towards throat; 15-18 subdigital lamellae on the 4th toe; and 21-23 ventral scales at midbody. The other, Cnemaspis thayawthadangyi sp. nov. is known only from the island group of Thayawthadangyi, in the Myeik Archipelago and is distinguished from all congeners by the possession of keeled pectoral scales; keeled ventral scales and abdominal scales; three precloacal pores; four femoral pores on each leg; smooth subcaudal scales, scales on median row enlarged; coloration of the gular region silver with dark-gray irregularly shaped streaks; 16-18 subdigital lamellae on the 4th toe; and 18-20 ventral scales at midbody. In addition, we address the taxonomic status of populations referred to as C. kandiana (a species now restricted to Sri Lanka) found in peninsular Thailand and provisionally assign them as C. cf. tanintharyi sp. nov. until further analysis can be conducted. Finally, we briefly discuss the biogeography of the South Asian clade of Cnemaspis.

Geographic structure of genetic variation in the Parachute Gecko Ptychozoon lionotum Annandale, 1905 across Indochina and Sundaland with descriptions of three new species.

An integrative taxonomic analysis of the Ptychozoon lionotum group across its range in Indochina and Sundaland recovers P. lionotum sensu lato Annandale, 1905 as paraphyletic with respect to P. popaense Grismer, Wood, Thura, Grismer, Brown, Stuart, 2018a and composed of four allopatric, genetically divergent, ND2 mitochondrial lineages. Multivariate and univariate analyses of continuous and discrete morphological and color pattern characters statistically and discretely diagnose each lineage from one another and together, with maximum likelihood and Bayesian inference analyses, provide the foundation for the recognition of each lineage as a new species-hypotheses corroborated with a Generalized Mixed Yule Coalescent species delimitation analysis. Ptychozoon cicakterbang sp. nov. ranges throughout Peninsular Malaysia to Pulau Natuna Besar, Indonesia; P. kabkaebin sp. nov. is endemic to northern and central Laos; and P. tokehos sp. nov. ranges from southern Thailand south of the Isthmus of Kra northward to Chiang Mai, fringing the Chao Phraya Basin and ranging southward through Cambodia to southern Vietnam. Ptychozoon lionotum sensu stricto ranges from northwestern Laos through southern Myanmar to eastern India. The phylogeographic structure within each species varies considerably with P. lionotum s.s. showing no genetic divergence across its 1,100 km range compared to P. cicakterbang sp. nov. showing upwards of 8.2% sequence divergence between syntopic individuals. Significant phylogeographic structure exists within P. tokehos sp. nov. and increased sampling throughout Thailand may require additional taxonomic changes within this species.

Filling the BINs of life: Report of an amphibian and reptile survey of the Tanintharyi (Tenasserim) Region of Myanmar, with DNA barcode data.

Despite threats of species extinctions, taxonomic crises, and technological advances in genomics and natural history database informatics, we are still distant from cataloguing all of the species of life on earth. Amphibians and reptiles are no exceptions; in fact new species are described nearly every day and many species face possible extinction. The number of described species continues to climb as new areas of the world are explored and as species complexes are examined more thoroughly. The use of DNA barcoding provides a mechanism for rapidly estimating the number of species at a given site and has the potential to record all of the species of life on Earth. Though DNA barcoding has its caveats, it can be useful to estimate the number of species in a more systematic and efficient manner, to be followed in combination with more traditional, morphology-based identifications and species descriptions. Herein, we report the results of a voucher-based herpetological expedition to the Tanintharyi (Tenasserim) Region of Myanmar, enhanced with DNA barcode data. Our main surveys took place in the currently proposed Tanintharyi National Park. We combine our results with photographs and observational data from the Chaung-nauk-pyan forest reserve. Additionally, we provide the first checklist of amphibians and reptiles of the region, with species based on the literature and museum. Amphibians, anurans in particular, are one of the most poorly known groups of vertebrates in terms of taxonomy and the number of known species, particularly in Southeast Asia. Our rapid-assessment program combined with DNA barcoding and use of Barcode Index Numbers (BINs) of voucher specimens reveals the depth of taxonomic diversity in the southern Tanintharyi herpetofauna even though only a third of the potential amphibians and reptiles were seen. A total of 51 putative species (one caecilian, 25 frogs, 13 lizards, 10 snakes, and two turtles) were detected, several of which represent potentially undescribed species. Several of these species were detected by DNA barcode data alone. Furthermore, five species were recorded for the first time in Myanmar, two amphibians (Ichthyophis cf. kohtaoensis and Chalcorana eschatia) and three snakes (Ahaetulla mycterizans, Boiga dendrophila, and Boiga drapiezii).

Open access genomic resources for terrestrial arthropods.

Genome sequencing initiatives like the Arthropod i5k project and other biodiversity genomics research rely on access to high quality DNA and/or tissue. Global collection initiatives such as the Smithsonian Global Genome Initiative (GGI) and its partner network, the Global Genome Biodiversity Network (GGBN) aim to provide access to these resources at high-quality standards. Here, we review progress toward providing genomic resources (tissues, DNA, genomes) for terrestrial arthropods, a megadiverse animal group, and compare progress in genome sequencing to all other animals.

Troublesome Trimes: Potential cryptic speciation of the Trimeresurus (Popeia) popeiorum complex (Serpentes: Crotalidae) around the Isthmus of Kra (Myanmar and Thailand).

The taxonomic identity of the Trimeresurus (Popeia) popeiorum complex from the Isthmus of Kra and to the north was investigated. Several studies over the last decade have produced several specimens and associated mtDNA sequence data for a variety of individuals of the T. popeiorum and "T. sabahi" complexes. Here, we combine four mitochondrial genes (12S, 16S, ND4, and CytB) from all available specimens in GenBank with the addition of five new specimens collected from the mainland, Tanintharyi Region of Myanmar. Maximum Likelihood and Bayesian analyses identified that T. popeiorum sensu lato is paraphyletic with two geographically distinct clades: a northern clade representing populations from northern Myanmar, Laos and northern Thailand and a southern clade representing samples from the Tanintharyi Region and adjacent west Thailand. While the two clades have considerable genetic distance, they appear to be morphologically identical, leading to the hypothesis that the southern clade represents a cryptic, undescribed species. Because they appear to be cryptic species and the limitation of only five specimens from the southern lineage, this does not permit us to formally describe the new species. In accordance to past molecular studies, we uncovered paraphyly and lack of genetic support for the validity of taxa within the T. sabahi complex. However, we suggest recognizing these populations as subspecies within T. sabahi.

Data Release: DNA barcodes of plant species collected for the Global Genome Initiative for Gardens Program, National Museum of Natural History, Smithsonian Institution.

The Global Genome Initiative has sequenced and released 1961 DNA barcodes for genetic samples obtained as part of the Global Genome Initiative for Gardens Program. The dataset includes barcodes for 29 plant families and 309 genera that did not have sequences flagged as barcodes in GenBank and sequences from officially recognized barcoding genetic markers meet the data standard of the Consortium for the Barcode of Life. The genetic samples were deposited in the Smithsonian Institution's National Museum of Natural History Biorepository and their records were made public through the Global Genome Biodiversity Network's portal. The DNA barcodes are now available on GenBank.

How many species and under what names? Using DNA barcoding and GenBank data for west Central African amphibian conservation.

Development projects in west Central Africa are proceeding at an unprecedented rate, often with little concern for their effects on biodiversity. In an attempt to better understand potential impacts of a road development project on the anuran amphibian community, we conducted a biodiversity assessment employing multiple methodologies (visual encounter transects, auditory surveys, leaf litter plots and pitfall traps) to inventory species prior to construction of a new road within the buffer zone of Moukalaba-Doudou National Park, Gabon. Because of difficulties in morphological identification and taxonomic uncertainty of amphibian species observed in the area, we integrated a DNA barcoding analysis into the project to improve the overall quality and accuracy of the species inventory. Based on morphology alone, 48 species were recognized in the field and voucher specimens of each were collected. We used tissue samples from specimens collected at our field site, material available from amphibians collected in other parts of Gabon and the Republic of Congo to initiate a DNA barcode library for west Central African amphibians. We then compared our sequences with material in GenBank for the genera recorded at the study site to assist in identifications. The resulting COI and 16S barcode library allowed us to update the number of species documented at the study site to 28, thereby providing a more accurate assessment of diversity and distributions. We caution that because sequence data maintained in GenBank are often poorly curated by the original submitters and cannot be amended by third-parties, these data have limited utility for identification purposes. Nevertheless, the use of DNA barcoding is likely to benefit biodiversity inventories and long-term monitoring, particularly for taxa that can be difficult to identify based on morphology alone; likewise, inventory and monitoring programs can contribute invaluable data to the DNA barcode library and the taxonomy of complex groups. Our methods provide an example of how non-taxonomists and parataxonomists working in understudied parts of the world with limited geographic sampling and comparative morphological material can use DNA barcoding and publicly available sequence data (GenBank) to rapidly identify the number of species and assign tentative names to aid in urgent conservation management actions and contribute to taxonomic resolution.

Rapid forest clearing in a Myanmar proposed national park threatens two newly discovered species of geckos (Gekkonidae: Cyrtodactylus).

Myanmar's recent transition from military rule towards a more democratic government has largely ended decades of political and economic isolation. Although Myanmar remains heavily forested, increased development in recent years has been accompanied by exceptionally high rates of forest loss. In this study, we document the rapid progression of deforestation in and around the proposed Lenya National Park, which includes some of the largest remaining areas of lowland evergreen rainforest in mainland Southeast Asia. The globally unique forests in this area are rich in biodiversity and remain a critical stronghold for many threatened and endangered species, including large charismatic fauna such as tiger and Asian elephant. We also conducted a rapid assessment survey of the herpetofauna of the proposed national park, which resulted in the discovery of two new species of bent-toed geckos, genus Cyrtodactylus. We describe these new species, C. lenya sp. nov. and C. payarhtanensis sp. nov., which were found in association with karst (i.e., limestone) rock formations within mature lowland wet evergreen forest. The two species were discovered less than 35 km apart and are each known from only a single locality. Because of the isolated nature of the karst formations in the proposed Lenya National Park, these geckos likely have geographical ranges restricted to the proposed protected area and are threatened by approaching deforestation. Although lowland evergreen rainforest has vanished from most of continental Southeast Asia, Myanmar can still take decisive action to preserve one of the most biodiverse places on Earth.

Resurrection of Bronchocela burmana Blanford, 1878 for the Green Crested Lizard (Squamata, Agamidae) of southern Myanmar.

Recent fieldwork in southern Tanintharyi revealed the presence of a small Green Crested Lizard in the wet evergreen forest. We generated mtDNA sequence data (ND2) that demonstrates that this population's nearest relative is Bronchocela rayaensis Grismer et al., 2015 of Pulau Langkawi, northwestern Peninsular Malaysia and Phuket Island. Morphologically the Burmese Bronchocela shares many features with Bronchocela rayaensis, which potentially would make this recently described Thai-Malay species a synonym of Bronchocela burmana Blanford, 1878; however, we interpret the genetic and morphological differences to reflect evolutionary divergence and recommend the recognition of both species.